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Abstract 

We calculate the number of photons produced by the parametric resonance 

in a cavity with vibrating walls. We consider the case that the frequency 

of vibrating wall is nui(n = 1,2,3,...) which is a generalization of other 

works considering only 2u>i, where uj\ is the fundamental-mode frequency of 

the electromagnètic field in the cavity. For the calculation of time-evolution 

of quantum fields, we introduce a new method which is borrowed from the 

time-dependent perturbation theory of the usual quantum mechanics. This 

perturbation method makes it possible to calculate the photon number for 

any n and to observe clearly the effect of the parametric resonance. 
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I. INTRODUCTION 



Recently, the photon creation in an empty cavity with moving boundaries, so called 
the dynamical Casimir effect QX|] has attracted much attention especially in a vibrating 
cavity 10-0]. It was also proposed that the high-Q electromagnètic cavities may provide a 
possibility to detect the photons produced in the dynamical Casimir effect [|7|||. Therein, 
they considered the vibrating wall with the frequency íl = 1u)\ and found the resonance 
excitation of the electromagnètic modes. 

For the solutions of Mathieu differential equation 

x + /x 2 (l + ecosfit)x = 0, (1.1) 

it is well known that parametric resonances occur when the frequency íl with which the 
parameter oscillates is close to any value 2/i/n with n integral 0. Since there are many mode 
freqüències for the differential equation describing the electromagnètic fields in a cavity, it is 
natural to consider other freqüències with Q = tvjJx- Now the method introduced in Ref. || 
cannot be used in these general cases because the generating function can be made only for 
the special case n — 2. Therefore we introduce a new perturbation method which leads to 
easily solvable equations. In principle, we may solve the equations to any order of expansion 
parameter e, but we calculate up to the first order of e, which is sufficient to see the effect 
of the parametric resonance. 

The organization of this paper is as follows. In Sec. || we review the scheme of the 



field quantization in the case of moving boundaries. In Sec. f I| we introduce the new 
perturbation method to find time evolution of quantum electromagnètic field. Here we 
write the dominant part of the solution of wave equation which results from the parametric 
resonance. In Sec. [IV] we calculate the number of photons created by the vibration of the 
boundary. The last section is devoted to the summary and discussion. Here, we discuss 
the physical properties of the parametric resonance in the coupled differential equations, 
and we estimate the photon number created in realistic situation. Finally, the higher order 
calculations are considered briefly. 



II. QUANTUM ELECTROMAGNÈTIC FIELDS IN A CAVITY WITH MOVING 

WALLS 



Let us consider an empty cavity formed by two perfect conducting walls, one being at 
rest at x = and the other moving according to a given law of motion L(t) when < t < T. 
The field operator in the Heisenberg representation A(x, t) associated with a vector potential 
obeys the wave equation (c = 1) 

d 2 A d 2 A 



dt 2 dx 2 

and can be written as 



(2.1) 



A(x,t) = J2KMx,t) + bl i r n (x,t)}. (2.2) 

n 

Here b\ and b n are the creation and the annihilation operators and ip n (x,t) is the corre- 
sponding mode function which satisfies the boundary condition if) n (0, t) = = if) n (L(t),t). 

For an arbitrary moment of time, following the approach of Refs. [|ÏT)| -|Ï2"|| , we expand the 
mode function as 

ip n (x, t) = Y^ Qnk(t)ip k (x, t) (2.3) 
k 

with the instantaneous basis 



»(*,£(»)) = ^ "i" (2-4) 



Here Q n k(t) obeys an infinite set of coupled differential equations 



Qnk + UJ k'(t)Qnk ~ 2A ^ gkjQnj + A ^ QkjQnj 

j j 



where X = L/L and 



+^ 2 J29jk9jiQni (2.5) 

3,1 



9kj = < ■ (2.6) 



(j = k) 



and the time-dependent mode frequency is 

«.(*) = gj- (2-7) 

For t < 0, the right hand side of Eq. ( |2.5|) vanishes and the solution in this region is 
chosen to be 

Then the quantum field can be written as 

í < 0) = Y,[b n -z==<p n & L o) + H.c] (2.9) 

n V 

L is the initial distance between the walls and uo n = j^. From the form of the Hamiltonian 
H = J2 n ^n(bnb n + |), we can intèrpret ò^ò n as the number operator associated with the 
photon with the frequency u n . 

After the oscillation of the wall (t >T), the solution of Eq. fl2.5|) with the initial condition 
(|2.8|) , can be written as 

Qnk{t >T)= a nk -= + (3 nk ^==. (2.10) 

V ^k V 

From ( |2.2j ) and (J2.3[), we have 

t >T) = J2[a n ^==<fn(x, L ) + H.c], (2.11) 

n V ^„ 

where the new creation and annihilation operators and a n are written in terms of and 
6 n , using the following Bogoliubov transformations 



a k = Y^n^nk + bÍ(3* nk ] 
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YX h Í<k + b n (3 nk \. (2.12) 



Further, it follows from H = Y, n UJ n( a n a n + è) ^hat the new number operator at 

t > T. 

If we start with a vacuum state |0&) such that b |0&) = 0, the expectation value of the 
new number operator is 



N k = (0 6 |4a fc |0 6 > = ^|^| 2 , (2.13) 

n=l 

which can be interpreted as the number of created photons. (One should note that the 
quantum state does not evolve in time in the Heisenberg picture.) 

III. TIME EVOLUTION OF QUANTUM ELECTROMAGNÈTIC FIELDS IN A 
CAVITY WITH AN OSCILLATING WALL 



In this section we find the time evolution of quantum field operator ( |2.2j ) by solving 
Eq. ( |2.5|) with the motion of the wall given by 

L(t) = Lq[1 + esin(m)]. (3.1) 

Here íl = jui —jtt/Lq and e is a small parameter characterized by the displacement of the 
wall. This is a generalization of the previous work in Ref. |§ where the special case 7 = 2 
was treated. For e < 1, having in mind that À(í) ~ e and taking the first order of e in the 
mode frequency (|2.7|) 

Uk (t) = — [1 + esin(fií)]" 1 , (3.2) 
Lo 

we can replace Eq. ( |2.5| ) by a pair of coupled first-order differential equations 

Qnk Pnk 

Pnk = -^fc 2 (l - 2e sin üt)Q nk + 2y ^ g k jP n j 

L i 

+TY,9k j Qnj + 0(e 2 ). (3.3) 

Let us now introduce 

X n , kT = ^ {Q nk ± i^*) (3.4) 

or inversely, 
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Then (fO) reads 



Qnk ^ 
Pnk ï 



2 



[Xn,k- + ^n,fc+] 



=F eVL COS fií ^ 9kj\l — h^nj- + -Xn 



i+J 



z 1 
=F -efi 2 sin fit 5- fei + X nJ+ ] 



Introducing the infinite dimensional column vector 



X n {t) 



X n ,i+ 

X n ,2- 



V • / 

the above equation can be written matrix form 

d 



dt 



X n (t) = V {0) X n (t) + eV {l) X n (t) 



and V<® and V® are matrices. The components of the matrices are 



and 



where 



v k%a' = iu k a5 kj 5 



T/í 1 ) _ V (,ni) s P si 1^ 
s=± 



3 1*', 7 \ fej. 



with s, a, o' = +, — . Here we used Q = jui and Uf. = kuj\ 



To find the solution of Eq. ( |3.8|) , we introduce a perturbation expansion: 

X n = X^ + eX^ + e 2 X^ + ---. (3.12) 

The iteration method used to solve this problem is similar to what we did in time-dependent 
perturbation theory. By inserting ( |3.12|) into Eq. (|3.8|) , identifying powers of e yields a series 
of equations: 

^CO) = y (0)^(0) (3 . 13) 

±XV> = VWX!® + V®Xg-\ (3.14) 
dt 



From the initial condition ( |2.8|) , we have the solution to zeroth order equation (|3.13| ) 

X ( X = tnk^-e-^ 1 . (3.15) 
Further, the Eq. (|3.14 ) is easily solved 



*SL(f) = ^ j dt'e-^' J2 VQ^X^. (3.16) 
Using ( |3.15| ) and ( |3.1Ü| ), it can be written explicitly as 

° J>' 

= üJxe™^ 1 * i' dt\v^ n _e- i{ak+ ^ +n)uJlt ' 
Jo 

+^ CT n-e +í( " <Tfe+7 "" ) " lí '), (3.17) 



that is, 



X nj k+(t) — v k _^ n _E n + ry+k (t) v k _ n _E n fe 7+fc (t), 



X^_(í) = v k _ tn _E* +J _ k (t) + v£_ tn _E%_^_ k (t), (3.18) 



where 



cuiíe lfca;ií , for m = 0, 

= < (3-19) 



Therefore we have, using Q3.18 ) and (|3] 



Q nk (t) = -L·e'^Xk 

V ZuJ k 

+ v k _^_E k n+1 _ k {t) + % + _, n _^_ 7 _ fc (í)] 

+ 0(e 2 ). (3.20) 

One should note that Q^ k includes terms proportional to iú\t which are the effects of para- 
metric resonance. In the usual situation, since uj\t ^> 1, only the resonance terms are 
dominant and the solution (|3.20|) becomes by retaining only them: 

Q nk (t) « -±=e-^& nk 
\/ lu k 

+ v k _ >n _e-^8 k , n+1 + 

(3.21) 



IV. NUMBER OF PHOTONS CREATED BY THE PARAMETRIC RESONANCE 



After some time interval T the wall stops at x = L Q . Then the wave function becomes 



Oínk 



\f%J k r y/^k 



ip k {x). 



(4.1) 



The mode function ( |2T3"| ) and its time-derivative should be continuous at t = T : 



Qnk(T)(pk = a nk—f== + Pnk 

k k \ V lUJ k 

Y \Qnk(T)(pk + Qnk(T)ípf 



D iüj fe T 



Wk 



2^ -iu) k a nk —== 

k V V ZuJ k 



+ iu k /3 nk 



(4.2) 



where we used ( |2.1(J| ) for the mode function at t > T. By multiplying <pi to both sides of the 
above equations and integrating, we can get a n k and /3 nk , which are 



S 



Olnk = ïUkQnk — Qnk + ~T X! 9klQnl 



L i ) W^k 

ftnk = \ ÍUkQnk + Qnk ~ ~jr dklQnlj ~pj^ z · (^-3) 

Retaining only the dominant terms (u)\T ^> 1) 

(3 nk = -e^Tvt^Jk^-ne-^. (4.4) 
Using ( |2.6| ) and (|3.11|) , finally we have 



IPnkf = ^nk{euiT) 2 8 ka - n . (4.5) 
Therefore the total number of photons created in the k th mode from the empty cavity is 

iV fc = ]T |/? nfe | = ^ (4.6) 
n=1 i 0, otherwise. 

This result is a generalization of Ref. || in the short time limit {euj\T <C 1) and it agrees 
with that result for 7 = 2 and k — 1. It should also be noted that the maximal number of 
photons are created at the mode frequency 

k = 1 or u k = ^. (4.7) 

for 7 = even and at its nearest neighbor freqüències k — (7 ± l)/2 for 7 = odd. 



V. DISCUSSION 



We changed the second order coupled differential equation (|2.5|) to the first order differ 



ential equation Q3.6D by introducing the new variables ( |3.4j ). This makes it easy to deal with 
the differential equation and to find the perturbation series by virtue of diagonalization of 
V^°\ The e 1 -order solution is found explicitly and it includes the terms proportional to time 
and this term is relatively large compared to other terms which include only oscillating parts. 
Considering only those dominant terms, we calculated the number of photons created after 
stopping of the wall vibration. The results show that the effect of parametric resonance is 



the largest at the half of the frequency of the vibrating wall (u>k = íl/2). This can be under- 
stood by considering the Mathieu equation ( |1.1| ) where the parametric resonance takes place 
most strongly for íl = 2/x. While, in the case 7 > 2, we see the other resonance effects in 
addition to Uk — 0/2, which is due to the effect of couplings with other mode freqüències in 
the cavity. This can be interpreted as the parametric resonance for the coupled differential 
equation. 

Following the discussions of Ref. [§], we estimate the rate of photon generation in the 
several modes. Our results (|4.6f) which are vàlid only in the limit euj\T <C 1 may be used 
to estimate the relative photon numbers depending on the mode frequency. For a long 
time, we assume that the photon numbers are proportional to time 0,0] or to exponential 
function of time [|]]. Then we have a photon distribution proportional to Eq. ( [4.6|) . For 
example, taking 7 = 4 we have N 2 = 4ÍVd and Ni = N 3 = 3Nn where is the number of 
created photons with Uf.- Here we used the same experimental parameters as considered in 
Ref. 0: e max ~ v s ó max /j7ic ~lx 10~ 8 , v s 6 max ~ 50 m/s, Ui/2ti ~ 10 GHz for L ~ 2 cm, 
and ÍVd ~ 300 photons during T = ls. For these experimental parameters, we expect that 
N 2 = 1200 photons and Aq = N 3 = 900 photons. 

The higher order calculations can be performed by taking the differential equation to 
higher order in e and iterating the integral. Here we note that when we make the e p -order 
differential equation like ( p.6|) from ( |2.5|) we consider the higher order coming from ( |3.2|) . 
However if we limit the problem to the dominant terms which result from the parametric 
resonance, the problem is simple. Let us explain these by considering the e 2 -order. The 
differential equation can be written, in the e 2 -order, as 



By integrating this differential equation it is clear that the term gives at most the term 
proportional to e 2 UiT. This is very small compared to {euj\T) 2 which comes from integrating 
the resonance term of Therefore it is enough to write 



(5.1) 



d 
dt 



(5.2) 
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to consider p-th order calculation. By iterations, we have 



pi 



0- p _l,...,(Tl i Sp,...,Sl 

W Sp X 
fc,cr;(Tp_i(Sp_i7-n),o-p_i 



s„_ i 

V X 
CTp-i(£p_i7-ra),CTp_i;<Tj,_2(£ í >_27-™)> CT í>-2 



X 



o"2 (S27— n) ,<T2 ;<ri (Ei 7— n) ,<ri 



(Ji (S17— n),(7i ;n,— ' 

(5.3) 

where = si + ... + Sfc. For large t, we should be careful about the convergency of the series. 
We expect that the convergency property may be provided by the factor l/p\. However the 
explícit calculations are too difficult to find the simple general form of the solution and this 



will be our concern of future work 13 
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